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ABSTRACT 

Chaichitehrani,  N.;  D’Sa,  E.J.;  Ko,  D.S.;  Walker,  N.D.;  Osburn,  C.L.,  and  Chen,  R.F.,  2014.  Colored  dissolved  organic 
matter  dynamics  in  the  Northern  Gulf  of  Mexico  from  ocean  color  and  numerical  model  results.  Journal  of  Coastal 
Research,  30(4),  800-814.  Coconut  Creek  (Florida),  ISSN  0749-0208. 

Colored  dissolved  organic  matter  (CDOM)  absorption  and  salinity  relationships  were  assessed  and  used  in  conjunction 
with  the  salinity  and  current  outputs  of  a  numerical  model  (Navy  Coastal  Ocean  Model  [NCOM])  to  study  CDOM 
dynamics  in  the  northern  Gulf  of  Mexico.  In  situ  CDOM  absorption  and  salinity  obtained  from  multiple  field  campaigns 
were  inversely  correlated  seasonally  (winter-spring  and  summer)  and  latitudinally  (inner-  and  outer-shelf  zones), 
suggesting  conservative  behavior  of  CDOM  distribution.  A  weaker  correlation,  during  summer  in  the  outer-shelf  zone, 
however,  indicated  stronger  effects  of  photooxidation  and  lower  masking  effects  from  riverine  CDOM.  Applying  these 
relationships  to  simulated  salinity  resulted  in  hourly  maps  of  CDOM  that  revealed  similarities  to  CDOM  patterns 
derived  from  SeaWiFS  satellite  imagery.  Further,  matchup  comparisons  between  model-derived  and  in  situ  CDOM 
absorption  were  statistically  sound  for  the  summer  (bias  =  —0.016,  root  mean  square  error  =  0.059,  r2  =  0.51  SI  =  0.28)  and 
the  winter-spring  periods  (bias  =  0.033,  root  mean  square  error  =  0.099,  r2  =  0.52,  SI  =  0.21).  Overlaying  the  model- 
derived  CDOM  maps  on  the  simulated  currents  revealed  the  strong  influence  of  currents  on  CDOM  advection. 
Downcoast  currents  during  the  nonsummer  months  led  to  persistent  advection  of  CDOM  westward  interrupted  by 
frequent  cold  front  events  that  flush  CDOM-laden  waters  out  of  the  coastal  bays  onto  the  inner  and  outer  continental 
shelves.  In  contrast,  the  upcoast  current  regime,  though  less  well  organized,  produces  a  more  significant  seaward 
advection  of  CDOM,  likely  due  to  the  Ekman  transport  and  subsequent  entrainment  by  mesoscale  eddies  over  the 
continental  slope. 


ADDITIONAL  INDEX  WORDS:  CDOM,  salinity,  SeaWiFS,  NCOM  model. 


INTRODUCTION 

Chromophoric  dissolved  organic  matter  (CDOM)  in  natural 
waters  absorbs  light  in  the  ultraviolet  (UV)  and  visible 
wavelength  range  and  plays  a  key  role  in  regulating  light 
attenuation  in  the  aquatic  medium.  In  coastal  waters,  CDOM  is 
predominantly  a  mixture  of  terrestrial  sources  ( e.g .,  humic  and 
fulvic  acids)  that  is  introduced  to  coastal  systems  by  river 
runoff.  By  contrast,  CDOM  in  the  open  ocean  is  predominantly 
a  mixture  of  autochthonous  sources,  such  as  phytoplankton 
and  its  subsequent  heterotrophic  processing  (Blough  and  Del 
Vecchio,  2002;  Nelson  et  al.,  2007,  and  references  therein). 
There  is  increasing  interest  in  monitoring  CDOM  due  to  its 
numerous  influences  and  applications.  Its  light-absorbing 
properties  can  significantly  impact  primary  productivity  by 
limiting  light  while  also  blocking  harmful  UV  radiation 
(Stedmon  et  al.,  2000).  CDOM  has  been  used  as  an  interme¬ 
diary  for  determination  of  dissolved  organic  carbon  concentra- 
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tions  using  remote  sensing  (Del  Castillo  and  Miller,  2008; 
Griffin  et  al.,  2011).  In  coastal  waters,  CDOM  has  been 
employed  as  a  semiconservative  water  mass  tracer,  as  an 
indicator  of  water  mass  mixing,  and  as  an  index  to  measure 
water  quality  (Blough  and  Del  Vecchio,  2002;  Coble  et  al.,  2004; 
D’Sa  and  Korobkin,  2008;  Ferrari  and  Dowell,  1998;  Le  Fouest 
et  al.,  2006).  Further,  CDOM  can  also  be  used  as  a  proxy  for 
retrieving  salinity  using  ocean  color  sensors  and  providing 
synoptic  maps  of  salinity  (Bowers  and  Brett,  2008;  D’Sa  et  al., 
2002).  Since  CDOM  plays  an  important  role  in  various 
biogeochemical  cycles  and  microbial  ecology  in  natural  waters 
(Moran  and  Zepp,  1997),  it  is  critical  to  study  its  distribution 
and  fate  in  the  aquatic  environment.  Although  satellite  ocean 
color  data  have  been  used  to  retrieve  and  monitor  CDOM  (D’Sa, 
Miller,  and  Del  Castillo,  2006;  Kahru  and  Mitchell,  2001; 
Mannino,  Russ,  and  Hooker,  2008),  their  more  regular  use  in 
monitoring  and  applications  is  hampered  by  cloud  coverage 
and  sun  glint  effects.  Cloud  coverage  often  limits  acquisition  of 
clear  imagery,  while  sun  glint  can  confound  precise  retrieval  of 
water’s  constituents  from  remote  sensors  (Yoder,  1999). 
Further,  use  of  satellite  data  is  limited  by  the  satellite’s 
temporal  and  spatial  resolution  and  the  time  of  the  satellite 
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pass.  In  contrast,  numerical  circulation  model  outputs  such  as 
temperature  and  salinity  that  can  be  directly  or  indirectly 
derived  from  satellite  data  can  be  generated  at  much  higher 
temporal  and  spatial  resolution.  Hence,  it  would  be  desirable  to 
find  an  approach  to  extract  CDOM  from  the  outputs  of  a 
numerical  model.  The  strong  linkage  between  salinity  and 
CDOM  reported  in  many  studies  especially  in  coastal  waters 
suggests  that  salinity  could  be  used  as  a  proxy  to  derive  CDOM 
from  simulated  salinity.  However,  in  order  to  use  numerical 
models  to  derive  an  ocean  color  product  such  as  CDOM  there  is 
a  requirement  for  (1)  a  strong  inverse  linear  relationship 
between  CDOM  and  salinity  and  (2)  a  reliable  numerical  model 
to  simulate  salinity. 

The  relationship  between  absorption  coefficients  of  CDOM 
(«cdom)  and  salinity  has  been  studied  for  different  regions 
(D’Sa  et  al.,  2002;  Ferrari  and  Dowell,  1998;  Stedmon  and 
Markager,  2003).  In  coastal  regions  highly  influenced  by 
CDOM-laden  river  runoff,  ccdom  often  declines  linearly  with 
increasing  salinity,  indicating  that  CDOM  behaves  conser¬ 
vatively  with  respect  to  salinity,  suggesting  a  simple  dilution 
process  with  neither  addition  nor  removal  of  terrestrial 
CDOM  (Blough  and  Del  Vecchio,  2002;  Blough,  Zafiriou,  and 
Bonilla,  1993;  Ferrari  and  Dowell,  1998;  Kowalczuk,  1999; 
Vodacek  et  al.,  1997).  Nonlinear  relationships  between 
acDOM  and  salinity  have  been  reported  in  different  aquatic 
environments,  including  coastal  regions  (Blough,  Zafiriou, 
and  Bonilla,  1993;  Chen  and  Gardner,  2004;  Gardner,  Chen, 
and  Berry,  2005;  Kowalczuk,  1999).  An  upward  or  a 
downward  deviation  from  an  inverse  linear  relationship 
between  CDOM  and  salinity  suggests  an  addition  (e.g., 
autochthonous  production)  or  removal  (e.g.,  photodegrada¬ 
tion,  adsorption)  of  CDOM  (D’Sa  and  DiMarco,  2009;  Uher  et 
al.,  2001)  or  the  mixing  of  multiple  water  masses  with 
different  CDOM  end  members  (Chen  and  Gardner,  2004; 
Ferrari  and  Dowell,  1998;  Stedmon  and  Markager,  2003; 
Uher  et  al.,  2001).  In  such  cases,  caution  should  be  taken  in 
using  salinity  as  a  proxy  for  CDOM. 

Currents  and  salinity  outputs  of  numerical  models  such  as 
the  Navy  Coastal  Ocean  Model  (NCOM)  can  be  used  in 
conjunction  with  CDOM  derived  using  CDOM-salinity  rela¬ 
tionships  to  investigate  physical  influences  on  CDOM  distri¬ 
butions.  Advection  of  CDOM  by  currents  is  an  important 
physical  process  that  could  affect  CDOM  distribution  in  the 
northern  Gulf  of  Mexico.  Jolliff  et  al.  (2003)  combined  a 
circulation  model  with  CDOM  photolysis  to  study  CDOM 
transport  over  the  west  Florida  shelf  (WFS).  Their  circulation 
model  included  advection-dispersion  equations  of  pollutant 
transport  with  modified  source  terms  to  include  CDOM 
transport.  Following  model  verification  with  in  situ  measure¬ 
ments  and  satellite  data,  simulation  results  were  employed  to 
study  CDOM  distribution  patterns  over  the  WFS.  Shulman  et 
al.  (2011)  combined  the  remotely  sensed  parameters  and  the 
outputs  of  the  NCOM  model  to  study  the  effects  of  the 
California  current  system  on  optical  parameters.  In  the 
northern  Gulf  of  Mexico,  better  insights  on  the  spatial  and 
temporal  distribution  and  fate  of  CDOM  could  be  obtained  by 
understanding  the  current  patterns  and  use  of  model  simula¬ 
tions. 


Various  factors  such  as  winds,  river  discharges,  and 
detached  eddies  from  the  Loop  Current  can  contribute  to 
the  generation  of  currents  over  the  Louisiana-Texas  shelf 
(Li,  Nowlin,  and  Reid,  1997;  Oey,  1995;  Walker,  2005;  Walker 
et  al.,  1996,  2005).  Results  from  a  modeling  study  indicate 
that  a  large  component  of  transport  over  the  inner  shelf  is 
produced  by  wind  (40-48%),  with  river  discharge  (28-33%) 
and  Loop  Current  (LC)  eddies  (19-33%)  accounting  for  the 
rest  (Oey,  1995).  Two  seasonal  current  patterns  that  depend 
on  wind  characteristics  are  distinguishable  over  the  shelf 
(Cho,  Reid,  and  Nowlin,  1998;  Cochrane  and  Kelly,  1986; 
Walker,  2005).  During  nonsummer  months  (from  September 
to  May),  the  large-scale  circulation  on  the  Louisiana-Texas 
continental  shelf  exhibits  a  cyclonic  gyre,  driven  mainly  by 
easterly  winds,  Ekman  transport,  and  geopotential  (Co¬ 
chrane  and  Kelly,  1986).  A  strong  coastal  jet  forms  the 
northern  limb  of  this  gyre,  which  is  modulated  by  discharge 
from  the  Mississippi  and  Atchafalaya  Rivers  (Cho,  Reid,  and 
Nowlin,  1998;  Cochrane  and  Kelly,  1986;  Walker,  2005)  in 
addition  to  winds.  Jarosz  and  Murray  (2005)  aptly  named  it 
the  Louisiana-Texas  Coastal  Current  (LTCC).  This  inner 
shelf  current,  also  known  as  a  downcoast  current,  is  well 
developed  with  a  clear  westward  direction  from  Louisiana  to 
Texas  and  southward  along  the  Texas  coast  (Cho,  Reid,  and 
Nowlin,  1998;  Walker,  2005;  Wiseman  and  Kelly,  1994).  The 
frequent  passage  of  cold  fronts  in  nonsummer  months 
subjects  the  LTCC  to  strong  W  and  NW  winds  for  one  to 
several  days  duration  (Walker,  2005).  The  wind  reversals 
have  a  rapid  impact  on  surface  circulation  and  often  on  the 
strength  and  duration  of  the  LTCC  (Murray  et  al.,  1998; 
Walker,  2005;  Walker  and  Hammack,  2000). 

Eastward-directed  outer  shelf  currents  form  the  outer  limb 
of  the  cyclonic  gyre  (Cochrane  and  Kelly,  1986)  and  the 
presence  of  LC  anticyclonic  eddies  often  enhance  the 
eastward  flow  on  the  outer  shelf  even  though  their  centers 
are  in  waters  ~1000  m  deep  (Cho,  Reid,  and  Nowlin,  1998; 
Walker,  2005;  Walker  et  al.,  1996).  The  downcoast  current 
weakens  and  disappears  during  summer  months  (June,  July, 
and  August)  (Cochrane  and  Kelly,  1986;  Walker,  2005)  as  a 
result  of  a  wind  shift  from  E  to  S  and  SW.  This  change  in  the 
wind  direction  starts  along  the  south  Texas  inner  shelf 
usually  in  May,  initiating  an  upcoast  current  that  strength¬ 
ens  and  expands  to  the  Louisiana  shelf  as  summer  progresses 
(Cochrane  and  Kelly,  1986;  Jarosz,  1997;  Walker,  2005).  The 
upcoast  current  is  not  as  well  developed  as  the  downcoast 
current  and  exhibits  a  disorganized  pattern  (Allahdadi,  Jose, 
and  Patin,  2012;  Crout,  Wiseman,  and  Chuang,  1984;  Jarosz 
and  Murray,  2005).  During  summer,  shelf-break  currents  are 
still  directed  to  the  east,  and  no  gyre  is  formed.  The 
downcoast  current  is  reestablished  by  September  when 
easterly  winds  resume  (Cochrane  and  Kelly,  1986;  Walker, 
2005),  and  it  persists  for  much  of  the  time  from  then  until 
sometime  in  May  (Cochrane  and  Kelly,  1986;  Walker,  2005). 
These  synoptic  and  seasonal  patterns  of  wind  behavior,  as 
well  as  the  seasonality  of  river  discharges  (Li,  Nowlin,  and 
Reid,  1997),  are  likely  to  contribute  strongly  to  the  formation 
of  distinct  patterns  of  CDOM  distribution  over  the  shelf. 

Therefore,  the  main  objective  of  the  present  study  is  to 
examine  physical  influences  on  surface  CDOM  distribution 
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Figure  1.  Map  of  the  study  region  in  the  northern  Gulf  of  Mexico 
encompassing  Mississippi,  Louisiana,  and  parts  of  Texas  coast.  Black  circles 
show  the  locations  for  which  MsLaTex  simulated  current  velocities  were 
acquired  for  March  and  June  2005.  MBD = Mississippi  Birdfoot  Delta,  SWP= 
Southwest  Pass,  BB  =  Barataria  Bay,  TB  =  Terrebonne  Bay,  AB  = 
Atchafalaya  Bay,  VB  =  Vermillion  Bay. 


using  NCOM  numerical  model  outputs  of  salinity  and  currents. 
The  study  includes  two  parts;  the  first  part  assesses  the 
relationship  between  CDOM  absorption  and  salinity  using  a 
large  database  of  in  situ  measurements  over  the  study  area  to 
derive  CDOM  from  simulated  salinity.  CDOM-salinity  rela¬ 
tionships  were  determined  for  different  zones  and  seasons, 
which  were  then  applied  to  the  simulated  salinity  maps.  The 
second  part  investigates  the  role  of  advective  influences  using 
simulated  current  patterns  on  CDOM  distribution  under  three 
main  physical  regimes  of  downcoast,  upcoast,  and  cold  front 
passages  in  the  northern  Gulf  of  Mexico. 

DATA  AND  METHODS 

Study  Area 

The  Northern  Gulf  of  Mexico,  comprising  both  the  inner  and 
outer  shelf  areas  between  88.2°W  to  95.5°W  and  latitude  27°N 
to  30.5°N,  were  studied  using  field  and  satellite  ocean  color  data 
along  with  the  outputs  from  a  numerical  model  (Figure  1). 

The  inner-shelf  region  of  this  area  is  characterized  by  high 
amounts  of  discharge  from  the  Mississippi  and  Atchafalaya 
Rivers,  which  is  characterized  by  distinct  seasonality  with  peak 
flow  in  spring  and  low  flow  in  autumn  (Muller-Karger  et  al., 
1991;  Walker  et  al.,  1994).  This  large  volume  discharge  of 


freshwater  strongly  affects  the  physical  and  biogeochemical 
processes  over  the  study  area  (D’Sa  and  Miller,  2003;  Lohrenz, 
Dagg,  and  Whitledge,  1990;  Muller-Karger  et  al.,  1991;  Walker 
and  Rabalais,  2006).  Although  the  river  discharge  increases 
CDOM  and  decreases  salinity  in  the  coastal  waters,  they  can 
both  be  modulated  by  shelf  hydrodynamics. 

Field  Data 

CDOM  and  salinity  data  were  obtained  from  nine  oceano¬ 
graphic  cruises  carried  out  in  2001  (April),  2005  (March,  May, 
July,  and  August),  2006  (January),  2007  (April  and  May),  and 
2009  (August)  (Table  1).  The  sampling  locations  are  shown  for 
the  winter-spring  (Figure  2 A)  and  the  summer  (Figure  2B) 
periods.  In  2005  (March,  May,  July,  and  August),  water 
samples  were  obtained  from  the  nearshore  coastal  waters 
strongly  influenced  by  the  Mississippi  and  Atchafalaya  Rivers 
(west  of  Southwest  Pass,  off  the  Terrebonne  Bay,  and 
Atchafalaya  shelf)  aboard  the  R/V  Gyre  (D’Sa  and  DiMarco, 
2009). 

Water  salinity  was  measured  using  a  conductivity-temper¬ 
ature-depth  profiler  (Sea-Bird  Electronic,  Inc.,  Bellevue, 
Washington,  U.S.A.)  during  all  cruises.  Water  samples  on  the 
R/V  Gyre  cruises  were  filtered  through  prerinsed  0.2- pm 
Nucleopore  membrane  filters.  Filtered  samples  were  kept 
refrigerated  at  4°C,  and  absorption  coefficients  were  measured 
on  single-beam  or  dual-beam  spectrophotometers.  Optical 
absorbance  spectra  measured  from  250  to  722  nm  using  Milli- 
Q  water  as  a  reference  were  further  corrected  for  baseline 
fluctuations  by  subtracting  the  average  absorbance  between 
715  and  722  nm  (D’Sa  and  DiMarco,  2009).  The  absorption 
coefficients  (m-1)  for  each  wavelength  were  then  obtained  using 
the  expression 

_  2.303A(i) 

aW  —  — i -  (i) 

where  A  is  the  absorbance,  1  is  the  wavelength,  and  l  is  the 
optical  path  length  in  meters  (0.1  m  or  0.5  m).  Similar  methods 
for  field  sampling  and  data  analysis  were  used  to  obtain  CDOM 
and  salinity  data  for  samples  acquired  from  other  field 
campaigns.  More  detailed  information  can  be  found  in  the 
Biological  and  Chemical  Oceanography  Data  Management 
Office  Web  site  (http://data.bco-dmo.org/jg/serv/BCO/ 
NACP_Coastal/GulfMexico/CDOM.html). 

Satellite  Data 

Data  from  the  SeaWiFS  ocean  color  sensor  (August  1997  to 
December  2010)  were  used  to  derive  CDOM  absorption 
coefficients  at  412  nm  (acd0m(412)).  Daily  SeaWiFS  Level- 1 
merged  local  area  coverage  data  with  1-km  resolution  were 
obtained  from  the  Ocean  Color  Web  site  run  by  the  National 
Aeronautics  and  Space  Administration  (NASA)  (oceancolor. 
gsfc.nasa.gov)  and  processed  using  SeaWiFS  Data  Analysis 
System  (SeaDAS)  into  Level-2  to  obtain  CDOM  absorption  by 
applying  the  (D’Sa  et  al.,  2006)  empirical  algorithm  given  by 
the  equation: 

«cdom(412)  =  0.227[i?(510)/J?(555)]-2'612  (2) 

where  Rrs(510)  and  Rrs(555)  are  remote  sensing  reflectances  at 
510  and  555  nm  wavelength  bands,  respectively.  The  robust- 
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Table  1.  Summary  of  cruises  undertaken  in  the  northern  Gulf  of  Mexico 
showing  periods  when  samples  were  acquired  for  this  study.  Average  values 
for  aCDOM(412)  and  salinity  shown  for  each  cruise. 


Cruise 

Date 

ocdom(412)  (m  1) 

Salinity  (psu) 

1 

12-16  April  2001 

0.23 

30.20 

2 

23-29  March  2005 

0.47 

27.07 

3 

20-25  May  2005 

0.43 

25.96 

4 

8-12  July  2005 

0.23 

29.27 

5 

18-24  August  2005 

0.25 

28.75 

6 

1-4  January  2006 

0.43 

23.50 

7 

7-10  May  2007 

0.87 

24.16 

8 

16-20  April  2007 

0.41 

31.75 

9 

18-20  August  2009 

0.34 

27.94 

ness  of  this  algorithm  was  demonstrated  based  on  matchup 
comparisons  with  field  data  (Chaichitehrani  et  al.,  2012, 
unpublished  data).  The  SeaWiFS- derived  surface  CDOM 
absorption  images  were  used  to  compare  with  NCOM-derived 
CDOM  maps. 

NCOM  Data 

A  high  spatial  resolution  (~1.9  km)  nested  Navy  Coastal 
Ocean  Model  (NCOM)  previously  developed  for  the  Mississip- 
pi-Louisiana-Texas  (MsLaTex)  coast  as  part  of  a  NASA- 
funded  project  was  used  in  this  study  (D’Sa  and  Ko,  2008). 
The  NCOM  MsLaTex  coastal  model  covers  the  region  between 
27°N  and  30.5°N  latitude  and  between  88.2°W  and  95.5°W 
longitude  (372  X  200  grid)  and  includes  the  Mississippi, 
Louisiana,  and  Texas  shelves  and  the  offshore  waters  in  the 
northern  Gulf  of  Mexico  (Figure  1).  The  MsLaTex  model  is 
nested  within  the  Intra-Americas  Sea  Nowcast/Forecast 
System  (IASNFS)  and  is  based  on  the  NCOM  (Barron  et  al ., 
2006;  Martin,  2000),  a  primitive  equation  model  with  hybrid 
sigma-z  vertical  coordinates.  The  model  horizontal  resolution  is 
between  1.8  and  1.9  km,  which  resolves  well  the  major 
circulation  features  on  the  shelf.  The  model  topography  is 
derived  from  NRL  2-min  global  topography  data  set  (DBDB2) 
and  improved  with  hydrographic  data  from  National  Geophys¬ 
ical  Data  Center.  There  are  111  rivers  or  freshwater  runoff 
points  included  in  the  model  with  daily  river  discharges  from 
the  Army  Corps  of  Engineers  and  U.S.  Geological  Survey 
stations  to  better  simulate  the  coastal  salinity.  The  surface 
forcing  of  wind,  heat,  and  air  pressure  are  from  the  Coupled 
Ocean/Atmosphere  Mesoscale  Prediction  System  high-resolu¬ 
tion  regional  weather  forecast  model.  The  boundary  conditions 
for  the  MsLaTex  model  were  imported  from  the  IASNFS. 
IASNFS  includes  the  Gulf  of  Mexico,  the  Caribbean,  and  parts 
of  the  western  North  Atlantic  with  horizontal  resolution  of  ~6 
km.  Both  IASNFS  and  MsLaTex  models  assimilate  sea  surface 
height  anomaly  data  from  satellite  altimeters  ( GFO ,  Jason-1, 
and  ERS-2)  and  sea  surface  temperature  data  from  the 
Advanced  Very  High  Resolution  Radiometer  to  reproduce  the 
important  mesoscale  Loop  Current  and  Loop  Current  Eddies  in 
the  Gulf  of  Mexico  (D’Sa  and  Ko,  2008;  Ko  et  al.,  2008). 


RESULTS 

CDOM-Salinity  Relationship 

Conservative  behavior  between  CDOM  absorption  and 
salinity  is  essential  to  employ  salinity  as  a  proxy  for  the 
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Figure  2.  Station  locations  of  cruises  conducted  during  (A)  winter-spring, 
and  (B)  summer  periods.  Zones  1  and  2  are  indicated. 


determination  of  CDOM.  For  the  study  area,  the  inverse 
linear  relationship  between  CDOM  absorption  and  salinity 
was  found  to  vary  seasonally  and  spatially;  seasonal  variabil¬ 
ity  in  discharges  from  the  Mississippi  and  Atchafalaya  Rivers, 
enhanced  water  column  mixing  during  the  winter-spring 
season,  and  elevated  CDOM  photobleaching  and  stratification 
during  summer  likely  contributed  to  this  variability.  Field 
data  collected  during  various  cruises  (Table  1)  were  exam¬ 
ined,  and  CDOM-salinity  relationships  were  established  for 
the  winter-spring  and  summer  periods  spatially  partitioned 
into  two  latitudinal  zones  (Figures  2A  and  B;  Table  2).  Zone  1 
comprises  the  shallow  areas  east  of  the  Birdfoot  Delta  and  the 
inner  shelf  west  of  the  delta.  Zone  2  comprises  the  deeper 
outer  shelf  extending  from  Mississippi  to  the  east  Texas  shelf. 
Zone  1  is  directly  affected  by  the  discharge  from  the 
Mississippi  and  Atchafalaya  Rivers.  In  zone  1,  effects  due  to 
photobleaching  are  most  likely  masked  by  the  lower  salinity 
waters  containing  elevated  levels  of  CDOM  (D’Sa  and 
DiMarco,  2009).  In  contrast,  zone  2  is  less  affected  by  the 
rivers’  freshwater,  and  therefore  less  conservative  behavior  of 
CDOM  and  salinity  is  expected.  This  is  generally  in 
accordance  with  the  shelf  division  suggested  by  Nowlin  et 
al.  (2005)  in  connection  with  studies  of  the  long-term  flow 
along  the  Texas-Louisiana  shelf. 

The  surface  CDOM  absorption  coefficient  at  412  nm  was  0.61 
±  0.57  m-1  at  an  average  surface  salinity  of  26.06  ±  8.22 
during  winter-spring  in  zone  1,  indicating  the  strong  river 
influence  on  CDOM  and  salinity  fields  in  the  nearshore  waters. 
During  the  same  period,  zone  2  exhibited  an  average  surface 
salinity  of  32.66  ±  4.47,  and  an  average  CDOM  absorption  of 
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Table  2.  Geographical  extension  of  zone  1  and  zone  2. 


Subset  Region 

Latitude  Range 

Longitude  Range 

Zone  1 

28.73°N-30.25°N 

88.00°W-94.00°W 

Zone  2 

28.00°N-28.729°N 

88.00°W-94.00°W 

0.203  ±  0.18  m_1,  which  was  three  times  less  than  that  of  the 
zone  1.  During  summer,  surface  CDOM  in  zone  1  (0.32  ±  0.22 
m-1)  was  almost  two  times  less  than  the  winter-spring  period, 
while  the  mean  surface  salinity  over  this  zone  increased 
slightly  to  27.47  ±  3.62.  Lower  CDOM  in  zone  1  during  the 
summer  could  be  attributed  to  the  low  flow  of  the  Mississippi 
and  the  Atchafalaya  Rivers  and  to  photobleaching  effects  aided 
by  stratification  and  high  solar  insolation  (Boss  et  al.,  2001; 
D’Sa  and  DiMarco,  2009).  Zone  2  during  summer  exhibited  the 
lowest  CDOM  absorption  (0.064  ±  0.103  m-1)  and  highest 
salinity  (32.84  ±  2.72),  suggesting  loss  of  CDOM  due  to 
photochemical  decomposition  (Chen  and  Gardner,  2004; 
Opsahl  and  Benner,  1998). 

For  each  season  and  each  zone,  the  CDOM  absorption 
coefficient  at  412  nm,  acDOM(412),  were  plotted  as  a  function 
of  salinity  (Figures  3A-D;  Table  3).  High  correlations  between 
CDOM  and  salinity  were  obtained  for  both  the  winter-spring 
and  summer  periods  in  zone  1  (r2  =  0.88,  n  =  117;  r2  =  0.92,  n  = 
39,  respectively)  (Figure  3A-C)  suggesting  conservative  mixing 
between  coastal  and  oceanic  waters  in  support  of  results  of  D’Sa 
and  DiMarco  (2009).  Similarly,  for  zone  2,  CDOM  and  salinity 
were  highly  correlated  during  the  winter-spring  period  (r2  = 


In  situ  salinity 


In  situ  salinity  (PSU) 


In  situ  salinity  (PSU) 


Figure  3.  Correlations  between  ucdom(412)  and  salinity  (A)  for  the  winter¬ 
spring  period  in  zone  1,  (B)  for  the  winter-spring  period  in  zone  2,  (C)  for 
summer  in  zone  1,  (D)  for  the  summer  in  zone  2. 


Table  3.  Results  showing  the  slope  (a),  intercept  (b)  and  correlation 
coefficient  (r2)  of  the  regression  equation a  between  8lcdom(412)  and  salinity 
for  the  two  zones. 


Subset  Region 

Season 

Number  of  Samples 

a 

b 

r2 

Zone  1 

Winter-spring 

117 

-0.066 

2.35 

0.88 

Zone  2 

Winter-spring 

22 

-0.039 

1.49 

0.93 

Zone  1 

Summer 

39 

-0.060 

1.98 

0.92 

Zone  2 

Summer 

25 

-0.028 

0.98 

0.54 

a  An  inverse  (negative)  linear  relationship  was  fitted  between  acDOM(412) 
(m-1)  and  salinity  in  the  form  of  acDOM(412)  =  a  X  salinity  +  b,  where  a  is 
the  slope  and  b  is  the  intercept. 


0.93,  n  =  22)  (Figure  3B),  likely  as  a  result  of  mixing  associated 
with  frontal  passages  during  the  period.  In  contrast,  during 
summer,  ucdom(412)  exhibited  low  correlations  with  salinity 
(Figure  3D),  and  many  of  the  data  points  fell  below  a 
conservative  mixing  line,  which  could  be  attributed  to 
photooxidation  and  loss  of  CDOM  during  the  summer. 

Spatial  Distribution  of  CDOM  Based  on  Modeled 
Salinity 

The  accuracy  of  NCOM  MsLaTex  model  for  simulating 
salinity  was  examined  by  comparing  modeled  and  in  situ 
salinity  for  available  salinity  data  (Figure  4).  Error  statistics 
such  as  bias,  root  mean  square  error  (RMSE),  scatter  index  (SI), 
and  correlation  coefficient  (r2)  were  calculated  using  the 
following  equations: 


1  N 

Bias=jyD(y;-*;) 


(3) 


RMSE  = 


N 


i  n 


(4) 


Figure  4.  Comparison  of  in  situ  salinity  with  simulated  salinity.  The 
diagonal  line  is  the  1-to-l  line. 
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where  yt  denotes  simulated  values,  xt  denotes  field  measured 
values,  y  is  the  average  of  simulated  value,  and  x  is  the 
average  of  in  situ  values.  N  is  number  of  samples.  Error 
statistics  (bias  =  0.28,  RMSE  =  3.56,  SI  =  0.12,  and  r2  =  0.73 
for  n  =  225)  indicated  relatively  satisfactory  performance  of 
the  model  in  simulating  salinity. 

CDOM-salinity  relationships  (Table  3)  were  applied  to 
modeled  salinity  to  obtain  model-derived  CDOM  absorption 
maps  C e.g .,  Figures  5A  and  B).  These  were  compared  with 
SeaW/FS-derived  CDOM  absorption  imagery  obtained  using 
the  D’Sa,  Miller,  and  Del  Castillo  (2006)  algorithm.  Based  on 
availability  of  clear  satellite  imagery  for  the  study  area, 
modeled  outputs  for  18  October  2005  and  23  July  2005  (1900 
h  UT)  were  selected  as  representatives  for  nonsummer  and 
summer  periods,  respectively.  Figures  5A-C  and  6A-C  show 
simulated  salinity  maps,  model-derived  CDOM  maps,  and 
SeaWiFS- derived  CDOM  absorption  for  both  the  selected 
times.  The  simulated  salinity  on  18  October  2005  (Figure  5A) 
shows  that  the  entire  study  area  except  some  regions  adjacent 
to  the  mouths  of  the  Mississippi  and  Atchafalaya  Rivers  were 
influenced  by  high  salinity  values  (~35  psu  on  average).  The 
low-salinity  regions  were  limited  to  areas  around  the  Mis¬ 
sissippi  Birdfoot  Delta,  Southwest  Pass,  Barataria,  Terre¬ 
bonne,  Atchafalaya,  and  Vermilion  Bays.  This  was  not 
unexpected,  since  river  discharges  are  normally  at  their  lowest 
in  autumn. 

The  inverse  and  linear  relationships  between  CDOM  and 
salinity  determined  for  the  study  area  implies  that  high  CDOM 
is  expected  in  the  low-salinity  coastal  region  and  low  CDOM  in 
the  high-salinity  offshore  waters.  This  is  illustrated  in  a  map  of 
surface  CDOM  absorption  derived  from  simulated  salinity 
(Figure  5B).  A  comparison  of  the  model-derived  surface  CDOM 
and  SeaWiFS- derived  CDOM  for  18  October  2005  at  1830 
hours  UT  (Figure  5C)  shows  similar  patterns  of  CDOM 
distributions.  Entrainment  of  CDOM  from  the  Southwest  Pass 
of  the  Mississippi  River  toward  the  inner  shelf  in  the  simulated 
map  appeared  similar  to  the  SeaW/FS-derived  image,  as  was 
the  pattern  at  the  mouth  of  the  Atchafalaya  Bay. 

Simulated  salinity  was  selected  for  23  July  2005  at  1900 
hours  UT  as  representative  of  summer  along  with  SeaWiFS 
imagery  for  the  same  day  (Figure  6).  The  simulated  salinity 
(Figure  6A)  depicts  low-salinity  filaments  extending  outward 
from  the  South  Pass  of  the  Mississippi  River  delta.  This  pattern 
appears  as  high-surface  CDOM  in  the  vicinity  of  South  Pass 
extending  toward  the  shelf  break  until  it  disperses  and  mixes 
with  high-salinity  offshore  waters.  Model  simulation  also 
shows  the  CDOM-laden  low-salinity  waters  spreading  east  of 
the  Birdfoot  Delta  toward  the  continental  shelf  of  the  NE  Gulf 
of  Mexico  (Figure  6B).  Although  the  derived  CDOM  exhibits 
high  consistency  with  the  satellite  image  (Figure  6C),  espe¬ 
cially  for  coastal  areas,  overestimates  in  satellite-derived 
CDOM  in  the  very  nearshore  waters,  especially  west  of  the 
Atchafalaya  Bay  system,  could  be  due  to  interference  from 
seawater  constituents  such  as  suspended  sediments  and  algal 
biomass  (D’Sa,  2008). 


The  ability  of  the  model  to  simulate  CDOM  was  investigated 
by  comparing  the  MsLaTex  model-derived  CDOM  with  in  situ 
measurements  for  the  winter-spring  and  the  summer  periods 
(Figures  7A  and  B,  respectively).  Model  data  used  for  matchup 
comparisons  were  acquired  from  simulations  for  May,  March, 
July,  and  August  2005  concurrent  with  in  situ  data  and  were 
independent  of  the  data  used  for  the  investigation  of  CDOM 
conservative  behavior.  Results  show  that  the  model-derived 
CDOM  was  quantitatively  in  agreement  with  the  field  data. 
For  the  winter-spring,  average  measured  and  modeled 
«cdom(412)  was  0.44  m-1  and  0.48  m-1  (bias  =  0.03;  RMSE  = 
0.099;  SI  =  0.21,  and  r2  =  0.52),  while  for  the  summer,  the 
average  values  of  ucdom(412)  were  0.2  m-1  and  0.18  m-1  (bias  = 
-0.016,  RMSE  =  0.059,  SI  =  0.28,  and  r2  =  0.51),  respectively. 
These  results  demonstrate  the  MsLaTex  model’s  ability  to 
predict  patterns  of  surface  CDOM  in  the  study  area.  In 
addition,  the  match-up  comparison  between  SeaWiFS- derived 
CDOM  using  the  D’Sa,  Miller,  and  Del  Castillo  (2006) 
algorithm  and  in  situ  measurements  for  the  winter-spring 
and  the  summer  periods  was  shown  by  Tehrani  et  at.  (2013). 

Effect  of  Advection  on  CDOM  Distribution 

Advection  by  currents  is  the  main  physical  process  that 
affects  the  distribution  of  CDOM  after  being  discharged  into 
the  marine  environment  by  rivers  (here  the  Mississippi  and  the 
Atchafalaya  Rivers)  (Boss  et  at .,  2001).  Seasonal  or  weather- 
related  reversals  in  the  direction  of  currents  on  the  shelf  along 
with  the  seasonal  variations  in  rivers’  discharge  and  CDOM 
levels  can  lead  to  significant  spatial  and  temporal  variabilities 
in  surface  CDOM  distribution.  In  order  to  study  CDOM 
distribution  in  relation  to  horizontal  advection,  surface  CDOM 
derived  from  modeled  salinity  was  superposed  over  the 
simulated  current.  This  revealed  that  the  CDOM  distribution 
patterns  corresponded  to  cold  front  passage  events  and 
expected  patterns  for  the  seasonal  downcoast  and  upcoast 
current  regimes  during  2005.  These  maps  suggest  strong 
physical  influences  on  CDOM  distribution  and  are  examined 
further  for  each  of  the  above  cases. 

Figure  8  depicts  hind-cast  CDOM  distributions  overlaid  with 
modeled  surface  currents  for  5  May  2005  (1800  h  UT)  that  show 
a  representative  downcoast  current  regime,  except  that  the 
offshore  limb  of  the  cyclonic  gyre  is  weak  to  nonexistent.  The 
general  current  direction  over  the  inner  shelf  is  westward,  with 
an  average  current  velocity  of  0.5  m/s.  Approaching  the  shelf 
break,  current  velocity  increases  to  approximately  1  m/s.  A 
large  anticyclonic  Loop  Current  Eddy  offshore  of  the  Mis¬ 
sissippi  Birdfoot  Delta  induced  strong  eastward  currents  at  the 
shelf  break.  In  general,  the  persistent  direction  of  shelf 
currents  to  the  west  entrained  CDOM  from  coastal  areas 
toward  the  inner  and  midshelf  regions.  Current  characteristics 
during  May  were  further  studied  by  examining  time  series 
NCOM  currents  at  points  1  and  2  (Figures  9A  and  D)  located 
south  of  the  Mississippi  Birdfoot  Delta  and  off  the  Terrebonne 
Bay,  respectively  (see  Figure  1  for  location).  At  point  1,  W  and 
SW  currents  were  dominant  for  the  first  5  days  of  May  (Figure 
9A),  while  SE  to  NE  currents  were  prevalent  for  the  rest  of  the 
month.  Point  2  exhibited  longer  periods  of  downcoast  current 
patterns  with  intermittent  southward  changes  in  current 
direction  induced  by  cold  front  passage  events.  CDOM-rich 
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Figure  5.  (A)  MsLaTex  simulated  salinity  map,  (B)  model-derived  CDOM  absorption  map  for  18  October  2005,  at  1900  h  UT,  and  (C)  SeaWfFS-derived  CDOM 
absorption  map  for  18  October  2005,  at  1830:24  h  UT.  (Color  for  this  figure  is  available  in  the  online  version  of  this  paper.) 


waters  (acDOM(412)  1.2  m-1  or  more)  west  of  the  Mississippi 
Birdfoot  Delta  and  off  the  Atchafalaya  Bay  were  advected  by 
westward  shelf  currents  downcoast.  The  high  surface  CDOM 
south  of  the  Birdfoot  Delta  was  not  entrained  toward  the  outer 
shelf,  since  the  current  pattern  implies  an  east  to  west  CDOM 


advection.  CDOM  absorption  was  relatively  high  (0.8-1  m-1)  in 
the  Louisiana  Bight  or  the  region  between  Southwest  Pass  of 
the  Mississippi  Birdfoot  Delta  and  Barataria  Bay.  The  modeled 
current  patterns  indicate  that  CDOM-laden  waters  were 
trapped  in  this  region  as  a  result  of  an  anticyclonic  (clockwise) 
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Figure  6.  (A)  MsLaTex  simulated  salinity  map,  (B)  model-derived  CDOM  absorption  map  for  23  July  2005,  at  1900  h  UT,  and  (C)  SeaWiFS- derived  CDOM  map 
for  23  July  2005,  at  1851:23  h  UT.  (Color  for  this  figure  is  available  in  the  online  version  of  this  paper.) 


gyre  generated  by  the  rotation  of  the  Mississippi  River  plume 
that  originated  off  Southwest  Pass  and  other  passes  to  the  east 
(Figure  10). 

Figure  11  illustrates  the  model-derived  CDOM  pattern 
overlaid  with  modeled  currents  on  11  March  2005  at  1200 


hours  UT,  associated  with  a  cold  front  passage  event.  Frequent 
cold  fronts  disturbed  the  westward  downcoast  current  system 
during  March  (Figures  9C  and  E).  For  example,  this  event  was 
associated  with  NW  winds  that  induced  S  to  SE  currents  with  a 
velocity  of  0.2  to  0.8  m/s  on  the  inner  Louisiana-Texas  shelf.  On 
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In  situ  acDOM(412)  (m_1) 


Figure  7.  Match-up  comparison  between  model-derived  CDOM  absorption 
and  in  situ  CDOM  (A)  for  the  winter-spring  period,  and  (B)  for  the  summer. 


the  outer  shelf,  the  currents  were  eastward,  most  likely  due  to 
the  impact  of  large  mesoscale  eddies.  Time  series  of  currents  for 
both  point  1  (Figure  9B)  and  point  2  (Figure  9E)  during  this 
period  showed  significant  southward  currents  as  a  result  of  cold 
front  outbreaks.  The  current  reached  a  maximum  velocity  of 
0.8  and  0.6  m/s  at  points  1  and  2,  respectively.  The 
corresponding  CDOM  pattern  was  strongly  tied  to  the  current 
pattern.  The  general  southward  direction  of  currents  on  the 
inner  shelf  advected  CDOM-laden  waters  toward  the  outer 
shelf.  CDOM-rich  waters  (CDOM  absorption  of  1.2  m-1  and 
greater)  from  the  Atchafalaya-Vermillion  and  Terrebonne- 
Timbalier  Bays  were  flushed  out  and  injected  onto  the  inner 
shelf  during  the  cold  front  passage  event.  The  same  process  also 
advected  CDOM  offshore  from  Texas  estuaries  and  coastal 
areas  onto  the  eastern  Texas  shelf. 


Typical  surface  currents  influencing  CDOM  distribution  for 
the  summer  (22  June  2005,  at  0000  h  UT)  are  illustrated  in 
Figure  12.  Shelf  currents  were  relatively  disorganized  with  a 
general  eastward  direction,  which  is  consistent  with  an  upcoast 
current  regime.  Shelf-break  currents  to  the  S  and  SW  of  the 
Birdfoot  Delta  were  eastward  and  influenced  by  a  mesoscale 
anticyclonic  eddy  located  offshore  of  the  Mississippi  Canyon  at 
latitude  27°30/N.  The  average  inner-shelf  current  velocity  at 
this  time  was  ~0.3  m/s.  Vector  representation  of  currents  at 
point  1  (Figure  9C)  during  June  showed  a  general  upcoast 
pattern  with  current  velocity  up  to  0.8  m/s.  There  was  a  general 
eastward  current  direction  for  the  first  9  days  of  June,  but  later 
strong  SW  currents  with  velocities  up  to  1.5  m/s  affected  the 
shelf  area  for  2  days.  The  direction  of  S  currents  changed  to  NE 
as  a  result  of  the  upcoast  current  regime  on  22  June.  As 
observed,  the  current  pattern  at  point  1  was  complex  and  could 
be  attributed  to  factors  such  as  wind  forcing,  Loop  Current 
eddies,  and  the  influence  of  the  Mississippi  River  plume 
(Martfnez-Lopez  and  Zavala-Hidalgo,  2009;  Nowlin  et  al ., 
2005;  Walker  et  al .,  2005).  Point  2  (Figure  9F)  experienced 
low  current  velocities  during  the  first  10  days  of  June  (less  than 
0.2  m/s);  within  the  next  7  days  strong  NW  to  NE  currents 
(maximum  velocity  of  1  m/s),  likely  induced  by  a  mesoscale 
eddy,  were  dominant  at  this  location.  The  current  velocities 
decreased  by  20  June  and  by  the  end  of  June  increased  and 
were  directed  southward  with  velocity  of  0.5  m/s  under  the 
influence  of  a  mesoscale  eddy  located  offshore  of  the  Mississippi 
Canyon  at  the  latitude  of  27o30rN.  On  the  eastern  Texas  shelf, 
elevated  CDOM  was  widely  distributed  over  the  shelf, 
presumably  as  a  result  of  Ekman  transport  induced  by  W  and 
SW  winds.  An  increase  in  current  velocity  (up  to  0.6  m/s)  off  the 
Atchafalaya  Bay  advected  CDOM-laden  waters  toward  the 
Louisiana  shelf,  while  an  increase  in  current  velocity  south  of 
Barataria  Bay  entrained  CDOM  southward. 

DISCUSSION 

CDOM  and  Salinity  Behavior 

CDOM  can  be  used  as  a  passive  proxy  to  track  freshwater 
masses  and  salinity  (Binding  and  Bowers,  2003;  D’Sa  and 
Miller,  2003;  D’Sa  et  al.,  2002;  Ferrari  and  Dowell,  1998).  The 
negative  linear  relationship,  detected  for  waters  within  our 
study  region,  between  salinity  and  CDOM  demonstrates  a 
reasonably  conservative  behavior  of  these  two  parameters 
modulated  by  seasonality  of  river  discharge.  The  dependency  of 
CDOM  on  salinity  was  demonstrated  for  two  different  zones 
during  the  winter-spring  and  the  summer  periods.  CDOM  and 
salinity  were  highly  correlated  in  zone  1  for  the  winter-spring 
and  summer  periods,  suggesting  CDOM  distribution  was 
regulated  by  mixing  (Chen  et  al.,  2007).  Elevated  levels  of 
CDOM  associated  with  low-salinity  waters  in  zone  1  indicate 
allochthonous  sources  of  CDOM  (Toming  et  al.,  2009).  In 
addition,  high  concentration  of  CDOM  introduced  by  the 
Mississippi  and  the  Atchafalaya  River  runoff  to  the  inner  shelf 
can  mask  variations  induced  by  removal  and  addition  processes 
such  as  photobleaching  and  microbial  activities,  respectively 
(Chen  et  al.,  2004;  D’Sa  and  DiMarco,  2009).  In  fact,  Walker  et 
al.  (2005)  suggested  that  as  much  as  75%  of  the  Mississippi 
River  discharge  may  be  advected  to  the  shelf  region  west  of  the 
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Figure  8.  Simulated  current  pattern  for  5  May  2005,  at  1800  h  UT  (color  bar  indicates  CDOM  absorption  coefficient  at  412  nm  (m  1)).  (Color  for  this  figure  is 
available  in  the  online  version  of  this  paper.) 


Mississippi  Birdfoot  Delta.  This  indicates  that  the  Terrebonne 
Bay  area  is  highly  affected  by  CDOM-rich  water  from  the 
Mississippi  River.  Solar  illumination  and  water  column 
stratification  during  the  summer  causes  photochemical  decom¬ 
position  of  CDOM,  which  disturbs  the  inverse  linear  relation¬ 
ship  between  CDOM  and  salinity.  However,  the  CDOM- 
salinity  relationship  for  zone  1  in  summer  exhibited  a  strong 
inverse  linear  relationship  ( i.e .  conservative  behavior).  This 
could  be  attributed  to  the  masking  effect  of  the  Mississippi  and 
the  Atchafalaya  discharge  (D’Sa  and  DiMarco,  2009).  Further¬ 
more,  the  required  time  scale  for  CDOM  photoxidation  is  days 
to  weeks,  which  is  longer  than  the  time  needed  for  river-laden 
CDOM  to  affect  the  inner  shelf  and  coastal  areas  (Hitchcock  et 
al.,  2004;  Vodacek  et  al.,  1997). 

A  different  seasonal  behavior  in  terms  of  the  CDOM-salinity 
relationship  was  observed  in  zone  2.  A  strong  conservative 
behavior  of  CDOM  and  salinity  during  the  winter-spring 
period  (r2  =  0.93)  was  observed.  In  contrast,  the  summer  time 
relationship  suggested  a  poor  conservative  behavior  of  CDOM 
(r2  =  0.54),  which  indicated  the  removal  of  CDOM,  presumably 
by  photodegradation.  The  footprint  of  CDOM  photodegradation 
can  be  identified  in  Figure  3D  by  data  points  that  were 
scattered  below  the  linear  CDOM-salinity  relationship.  The 
strength  of  photobleaching  over  this  offshore  zone  in  the 
summer  was  accentuated  by  the  low  river  discharge  and 
absence  of  its  masking  impact,  which  is  less  prominent  over 
zone  2  during  the  summer  time.  Dilution  of  river  CDOM 
increases  the  relative  depth  of  photodegradation,  enhancing 
CDOM’s  removal  in  surface  waters  (Osburn,  Retamal,  and 
Vincent,  2009).  In  addition,  an  increase  in  salinity  over  zone  2 
may  increase  CDOM’s  photoreactivity,  especially  at  longer 
wavelengths  used  for  remote  sensing  (Osburn  and  Morris, 
2003;  Osburn,  O’Sullivan,  and  Boyd,  2009). 

The  conservative  relationships  between  CDOM  and  salinity 
obtained  using  a  seasonal  and  zonal  approach  generally 
resulted  in  higher  accuracy  CDOM  concentration  maps  from 
simulated  salinity.  The  similarity  between  CDOM’s  distribu- 


(D) 


1-May  10-May  20-May  30-May 


(E) 


10- June  20-June  30-June 


Figure  9.  Time  series  of  current  (A)  at  point  1  for  May  2005,  (B)  at  point  1  for 
March  2005,  (C)  at  point  1  for  June  2005,  (D)  at  point  2  for  May  2005,  (E)  at 
point  2  for  March  2005,  and  (F)  at  point  2  for  June  2005  (see  Figure  1  for 
points  1  and  2  locations). 
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Figure  10.  Anticylonic  gyre  pattern  in  the  Lousiana  Bight  on  5  May  2005  (1800  h  UT)  (color  bar  indicates  CDOM  absorption  coefficient  at  412  nm  (m  1)).  (Color 
for  this  figure  is  available  in  the  online  version  of  this  paper.) 


tion  maps  obtained  from  model  and  SeaWiFS  data  suggests  that 
the  approach  for  producing  CDOM  synoptic  maps  using 
simulated  salinity  and  CDOM-salinity  conservation  concepts 
is  an  appropriate  method.  However,  discrepancies  between 
model-derived  and  SeaWiFS- derived  surface  CDOM  are  ob¬ 
served  in  estuaries  and  very  shallow  waters.  These  could  be  due 
to  errors  in  atmospheric  correction  and  interference  by  the  other 
optically  active  seawater  constituents  (D’Sa,  2008;  Le  Fouest  et 
al.,  2006;  Walker  and  Rabalais,  2006).  Model  deficiency  in 
simulating  salinity  in  coastal  areas  may  also  contribute  to  a 
lower  accuracy  of  CDOM  estimation  in  these  regions.  Model- 
derived  CDOM  maps  can  be  used  as  supplemental  data  to  fill 
temporal  and  spatial  gaps  in  the  satellite  data. 

Advection  and  CDOM  Distribution 

The  role  of  horizontal  advection  on  CDOM  distribution  on  the 
inner  and  outer  shelves  was  demonstrated  during  different 


current  regimes.  Cold  front  passages  induce  southward 
currents  and  strongly  contribute  to  the  CDOM  advection  over 
the  Louisiana  shelf.  The  current  can  flush  CDOM-laden  water 
out  of  the  bays  toward  the  shelf.  Walker  and  Hammack  (2000) 
estimated  that  30-50%  of  the  Atchafalaya  Bay  system’s  volume 
is  flushed  onto  the  inner  shelf  during  cold  front  passages.  Feng 
(2009)  reported  that  strong  cold  fronts  could  flush  more  than 
40%  of  the  bay  waters  to  the  shelf  within  a  period  of  less  than  40 
hours,  indicating  these  events  are  responsible  for  major 
outward  flushing  from  the  bays.  The  significant  effect  of  N 
and  NW  winds  on  the  flushing  of  water  out  of  the  bays  has  also 
been  demonstrated  using  satellite  data  and  numerical  models 
(Cobb,  Keen,  and  Walker,  2008;  Walker  and  Hammack,  2000). 
Thus,  relatively  high  southward  current  velocities  over  the 
shelf  associated  with  the  cold  fronts  also  advect  CDOM.  Figure 
13B  illustrates  CDOM  distribution  superposed  with  isobaths 
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Figure  12.  Current  pattern  for  22  June  2005,  at  0000  h  UT  superposed  on  CDOM  (color  bar  indicates  CDOM  absorption  coefficient  at  412  nm  (m  1)).  (Color  for 
this  figure  is  available  in  the  online  version  of  this  paper.) 


over  the  inner  and  outer  shelves  for  11  March  2005.  Owing  to 
the  influence  of  the  wind  and  southward  advection  caused  by 
the  cold  front,  elevated  CDOM  was  detected  seaward  to  the  60- 
m  isobath  off  the  Terrebonne  Bay,  to  the  30-m  isobath  off  the 
Atchafalaya  Bay,  and  to  the  40-m  isobath  off  the  eastern  Texas 
coast.  Westward-directed  downcoast  currents  contributed  less 
to  advection  of  CDOM  toward  the  outer  shelf  compared  with 
the  cold  front-induced  currents.  While  cold  front-induced 
currents  advected  CDOM  to  ~50-m  isobath  (on  average), 
downcoast  current  regimes  transport  CDOM  to  ~30-m  isobath 
(Figures  13A  and  B). 

The  persistent  westward/downcoast  current  regime  may  lead 
to  the  formation  of  a  relatively  high  CDOM  zone  parallel  to  the 
Louisiana  and  eastern  Texas  coast  due  to  freshwater  discharge 
from  rivers.  In  the  Louisiana  Bight,  a  clockwise  quasi¬ 
permanent  easterly  wind-driven  gyre  traps  freshwater  dis¬ 
charge  from  the  Mississippi  River  and  results  in  an  increase  in 
stratification  (Walker  and  Rabalais,  2006)  and  relatively  high 
CDOM  levels. 

Significant  increases  in  buoyancy  and  stratification  over 
the  inner  shelf  due  to  advection  of  the  Mississippi  River 
discharge  have  been  reported  by  Wiseman  et  al.  (1997). 
Although  the  upcoast  current  regime  is  not  as  organized  as 
that  of  the  downcoast  current,  CDOM  can  be  advected  farther 
offshore  during  upcoast  currents.  As  shown  in  Figure  13C, 
CDOM  is  entrained  seaward  to  the  60-m  isobath  off  the 
Terrebonne  Bay,  the  50-m  isobath  off  the  Atchafalaya  Bay, 
and  the  50-m  isobath  off  the  eastern  Texas  coast.  This 
distribution  could  be  due  to  Ekman  transport  induced  by  W 
and  SW  winds,  which  are  prevalent  during  the  upcoast 
current  regime.  In  addition,  mesoscale  eddies  can  influence 
currents  on  the  outer  shelf,  contributing  to  offshore  entrain¬ 
ment  of  CDOM.  Further,  during  summer,  the  general  pattern 
is  upcoast  currents  transporting  salty  water  from  the 
western  Gulf  of  Mexico  and  Texas  to  the  Louisiana  shelf 
(Li,  Nowlin,  and  Reid,  1997).  This  advection  of  salty  water 


may  also  disturb  the  conservative  relationship  between 
CDOM  and  salinity  in  summer. 

SUMMARY  AND  CONCLUSION 

The  combined  use  of  a  numerical  model  and  remote  sensing 
data  provides  a  valuable  tool  to  investigate  CDOM  dynamics 
with  high  spatial  and  temporal  resolution  in  the  northern  Gulf 
of  Mexico.  Accordingly,  an  approach  was  developed  to  construct 
synoptic  maps  of  surface  CDOM  absorption  for  the  region  using 
modeled  salinity.  The  main  requirement  for  implementing  this 
approach  was  converting  salinity  to  CDOM.  This  can  be 
successfully  accomplished  when  CDOM  covaries  inversely 
with  salinity.  Field  data  of  CDOM  and  salinity  from  different 
field  campaigns  were  used  to  examine  the  seasonal  and  spatial 
dependence  of  CDOM  on  salinity.  The  conservative  behavior 
was  met  for  both  the  inner-  and  outer-shelf  zones  during  the 
winter-spring  period,  while  the  summer  conservative  behavior 
could  only  be  met  over  the  inner- shelf  zone  due  to  the  masking 
influence  of  rivers’  runoff  on  the  destructive  effects  of 
photodegradation.  The  resulting  CDOM-salinity  relationships 
were  applied  to  simulated  salinity  to  produce  CDOM  maps  for 
different  physical  regimes  over  the  study  area.  The  NCOM- 
derived  surface  CDOM  distribution  maps  were  compared  with 
SeaWiFS- derived  CDOM  absorption  maps  (D’Sa,  Miller,  and 
Del  Castillo,  2006),  and  a  satisfactory  agreement  was  observed. 
Also,  the  remarkable  agreement  between  in  situ  and  NCOM- 
derived  acDOM(412)  as  shown  by  a  quantitative  evaluation 
reveals  the  model’s  ability  to  predict  CDOM  with  hourly 
temporal  resolution.  Some  disparities  that  were  seen  between 
these  two  maps,  especially  in  estuaries  and  shallow  water 
regions,  could  be  attributed  to  the  model’s  drawback  in 
simulating  salinity  in  these  regions.  Satellite  CDOM  estimates 
contaminated  by  other  optically  active  seawater  constituents 
could  also  induce  further  discrepancies.  Simulated  currents 
were  further  applied  to  investigate  CDOM  advection  over  the 
study  area  by  overlaying  simulated  currents  over  the  model- 
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Figure  13.  Surface  CDOM  absorption  map  derived  from  model  (A)  on  5  May  2005  at  1800  hours  UT  (B)  on  11  March  2005  at  1200  h  UT  and  (C)  on  22  June  2005  at 
0000  h  UT.  (Color  for  this  figure  is  available  in  the  online  version  of  this  paper.) 


derived  CDOM  maps.  The  impact  of  advection  on  CDOM 
distribution  was  demonstrated  through  a  detailed  comparison 
for  the  upcoast  and  the  downcoast  current  regimes  with  the 
corresponding  model-derived  CDOM  patterns.  Southward 
currents  induced  by  cold  front  passage  events  advect  CDOM- 


laden  waters  from  coastal  bays  to  the  inner  and  outer  shelf, 
while  the  persistent  westward  downcoast  current  regime 
contributed  less  to  CDOM  advection  toward  the  outer  shelf. 
The  downcoast  current  regime  produces  continuous  CDOM 
advection  within  the  LTCC  along  the  Louisiana  and  eastern 
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Texas  coastlines.  Seaward  CDOM  advection  into  deep  water 
during  an  upcoast  current  regime  can  be  significant  due  to 
Ekman  transport  and  the  influence  of  me  so  scale  eddies  on 
outer  shelf  currents. 

This  study  demonstrates  the  suitability  of  the  approach  to 
generate  high-resolution  CDOM  maps  based  on  simulated 
salinity.  Since  the  base  of  this  approach  is  the  conservative 
behavior  of  CDOM  and  salinity,  it  is  possible  to  reverse  it. 
Hence  salinity  can  be  obtained  using  satellite  data.  This  would 
be  desired  when  salinity  data  are  not  available  or  reliable. 
However,  the  study  area  is  a  highly  complex  and  dynamic 
environment,  and  uncertainties  associated  with  microbial  and 
photochemical  removal  mechanisms  and  biological  sources  of 
CDOM  need  to  be  incorporated  in  future  studies. 
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